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ABSTRACT 

Formamide (NH 2 CHO) has been proposed as a pre-biotic precursor with a key role in 
the emergence of life on Earth. While this molecule has been observed in space, most of 
its detections correspond to high-mass star-forming regions. Motivated by this lack of 
investigation in the low-mass regime, we searched for formamide, as well as isocyanic 
acid (HNCO), in 10 low- and intermediate-mass pre-stellar and protostellar objects. 
The present work is part of the IRAM Large Programme ASAI (Astrochemical Surveys 
At IRAM), which makes use of unbiased broadband spectral surveys at millimetre 
wavelengths. We detected HNCO in all the sources and NH 2 CHO in five of them. 
We derived their abundances and analysed them together with those reported in the 
literature for high-mass sources. For those sources with formamide detection, we found 
a tight and almost linear correlation between HNCO and NH 2 CHO abundances, with 
their ratio being roughly constant -between 3 and 10- across 6 orders of magnitude in 
luminosity. This suggests the two species are chemically related. The sources without 
formamide detection, which are also the coldest and devoid of hot corinos, fall well 
off the correlation, displaying a much larger amount of HNCO relative to NH 2 CHO. 
Our results suggest that, while HNCO can be formed in the gas phase during the cold 
stages of star formation, NH 2 CHO forms most efficiently on the mantles of dust grains 
at these temperatures, where it remains frozen until the temperature rises enough to 
sublimate the icy grain mantles. We propose hydrogenation of HNCO as a likely 
formation route leading to NH 2 CHO. 

Key words: astrochemistry - methods: observational - stars: formation - ISM: 
molecules - ISM: abundances 


1 INTRODUCTION 


* E-mail: ana®taurus.phys.s.u-tokyo.ac.jp 

t Based on observations carried out with the IRAM 30-m Tele¬ 
scope. IRAM is supported by INSU/CNRS (France), MPG (Ger¬ 
many) and IGN (Spain). 

t Based on analysis carried out with the CASSIS software. 


One of the major questions regarding the origin of life on 
Earth is whether the original chemical mechanism that led 
from simple molecules to life was connected to metabolism 
or to genetics, both intimately linked in living beings. For¬ 
mamide (NH2CHO) contains the four most important ele- 
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ments for biological systems: C, H, O, and N, and it has 
recently been proposed as a pre-biotic precursor of both 
metabolic and genetic material, sug gesting a common che m- 
ical origin for the two mechanisms JSaladino et al.ir2012h . 

Formamide was detected for the first time in space 
by iR.ubin et alJ (Il97ll 'l towards Sgr B2 and later in Orion 
KL. However, dedicated studies of NH 2 CHO in molecu¬ 
lar clouds have started only very recently, as its poten¬ 
tial as a key prebiotic molecule has become more evi¬ 
dent. These studies present observations o f formamide in 
a number of massive hot molecular cores llBisschop et al.l 
120071 . lAdande et al.l 1201 ill, the low-mass p rotostellar ob¬ 
ject IRAS 16293-2422 (iKahane et al.ll2013ll. and the out- 
flow shock region s L1157-B1 and B2 jMendoza et ^ I I 2 OI 4 
lYamaguchi et all I2OI2II . Its detection in comet Hale -Bopp 
has also been reported (iBockelee-Morvan et aLlbood l. For¬ 
mamide is therefore present in a variety of star-forming 
environments, as well as on a comet of the Solar System. 
Whether this implies an exogenous delivery onto a young 
Earth in the past is a suggestive possibility that needs more 
evidence to be claimed. 

Establishing the formation route(s) of formamide in 
space remains a challenge. Different chemical pathways have 
been proposed, both in the gas-phase (e.g. Redondo et al.l 
2014^ and on g rain surfaces (e.g. iRaunier et al.l |2004 
Jones et al.ll20lil l. The present work represents an effort to 


try to understand the dominant mechanisms that lead to 
the formation of formamide in the interstellar medium. In 
particular, it seeks to investigate the possible chemical con- 
ne ction between NH 2 CH O and HNCO, which was proposed 
bv iMendoza et ahl (l2014l l. To this aim, we have performed a 
homogeneous search of NH 2 CHO and HNCO in a represen¬ 
tative sample of 10 star-forming regions (SFRs) of low- to 
intermediate-mass type, since most of the formamide detec¬ 
tions so far reported concentrate on high-mass SFRs. This 
is the first systematic study conducted within the context of 
the IRAM Large Program ASAI (Astrochemical Surveys At 
IRAM; P.I.s: B. Lefloch, R. Bachiller), which is dedicated to 
millimetre astrochemical surveys of low-mass SFRs with the 
IRAM 30-m telescope. 

The source sample and the observations are described 
in Sects. [2] and O respectively. Section |4] presents the spec¬ 
tra and describes the analysis carried out to obtain the 
abundances of NH 2 CHO and HNCO. Sect. [5] compares the 
derived abundances with those found in the literature for 
other SFRs, and discusses the formation mechanisms that 
are favoured by our results. Our conclusions are summarised 
in Sect. [G] 


cated study of Complex Organic Molecules (COMs) in this 
source, including NH 2 CHO, was recently carried out by 
Ijaber et al.1 (I20l4 . 


3 OBSERVATIONS AND DATA REDUCTION 

The data presented in this work were acquired with the 
IRAM 30-m telescope near Pico Veleta (Spain) and con¬ 
sist of unbiased spectral surveys at millimetre wavelengths. 
These are part of the Large Programme ASAI, whose obser¬ 
vations and data reduction procedures will be presented in 
detail in an article by Lefloch & Bachiller (2015, in prep.). 
Briefly, we gathered the spectral data in several observing 
runs between 2011 and 2014 using the EMIR receivers at 
3 mm (80-116 GHz), 2 mm (129-173 GHz), and 1.3 mm 
(200-276 GHz). The main beam sizes for each molecular 
line analysed are listed in Tables IBII and IB2I The three 
bands were covered for most of the sources. For Cep E, ad¬ 
ditional observations were carried out at 0.9-mm (E330 re¬ 
ceiver), while just a few frequencies were covered at 2 mm. 
The Fourier Transform Spectrometer (FTS) units were con¬ 
nected to the receivers, providing a spectral resolution of 
195 kHz, except in the case of L1544, for which we used the 
FTS50 spectrometer, with a resolution of 50 kHz, to resolve 
the narrow lines (AV ~ 0.5 km s“^) that characterise this 
region. The observations were performed in wobbler switch¬ 
ing mode with a throw of 180”. 

The data were reduced with the package CLASS90 of 
the GILDAS software collection^ Through comparison of 
line intensities among different scans and between horizon¬ 
tal and vertical polarisations, the calibration uncertainties 
are estimated to be lower than 10 % at 3 mm and 20% in 
the higher frequency bands. After subtraction of the con¬ 
tinuum emission via hrst-order polynomial htting, a final 
spectrum was obtained for each source and frequency band 
after stitching the spectra from each scan and frequency set¬ 
ting. The intensity was converted from antenna temperature 
(Taut) to main beam temperature (TmM using the beam effi¬ 
ciencies provided at the IRAM web sitqj. In order to improve 
the signal-to-noise (S/N) ratio, the 2- and 1-mm ASAI data 
were smoothed to 0.5 km s“^, except in the case of L1544, 
for which we kept the original spectral resolution. 

For 116293, we used the TIMASSS spectral data ob¬ 
tained with the IRAM 30-m telescope at 1, 2, and 3 mm. A 
detailed description of the observations and a n overview of 
the dataset are reported in ICaux et al.l (I2OIIII . 


2 SOURCE SAMPLE 

Our source sample consists of 10 well-known pre-stellar and 
protostellar objects representing different masses and evolu¬ 
tionary states, thus providing a complete view of the var¬ 
ious types of objects encountered along the first phases 
of star formation. Their basic properties are fisted in Ta¬ 
ble [T] All of them belong to the ASAI source sample ex¬ 
cept one: the Glass 0 protobinary IRAS 16293-2422 (here¬ 
after 116293), whose millimetre spectral survey, TIMASSS 
(The IRAS16293-2422 Millimeter And Submillimeter Spec¬ 
tral Survey), was published by ICaux et al.l |20l3). A dedi¬ 


4 RESULTS 
4.1 Line spectra 

We searched for formamide (NH 2 GHO) and isocyanic acid 
(HNGO) in our dataset using the CASSIS softwar^f] and 
the Cologne Datab ase for Molecular Spectroscopy (CDM^; 
iMuller et al.ll200l] . 2005) to identify the fines. For NH 2 CHO, 

^ http://www.iram.fr/IRAMFR/GILDAS/ 

^ http://www.iram.es/IRAMES/mainWiki/Iram30mEfficiencies 
3 CASSIS has been developed by IRAP-UPS/CNRS 
(http://cassis.irap.omp.eu) 

^ http://www.astro.uni-koeln.de/cdms/ 
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Table 1. Source sample and their properties. 


Source 

R.A.(J2000) 

Dec.(J2000) 

Visr 

(kms“^) 

d 

(pc) 

M 

(M©) 

-^bol 

(Lo) 

Type* 

References 

ASAI 

LI 544 

05:04:17.21 

-1-25:10:42.8 

-1-7.3 

140 

2.7 

1.0 

PSC 

1,2,3 

TMCl 

04:41:41.90 

-1-25:41:27.1 

-1-6.0 

140 

21 

— 

PSC 

1,4 

Bl 

03:33:20.80 

-1-31:07:34.0 

-1-6.5 

200 

1.9 

1.9 

Class 0 

5,6 

L1527 

04:39:53.89 

-1-26:03:11.0 

-1-5.9 

140 

0.9 

1.9 

Class 0, WCCC 

1,7,8 

L1157-mm 

20:39:06.30 

-1-68:02:15.8 

-1-2.6 

325 

1.5 

4.7 

Class 0 

7,8 

IRAS 4A 

03:29:10.42 

-1-31:13:32.2 

-1-7.2 

235 

5.6 

9.1 

Class 0, HC 

7,8 

SVS13A 

03:29:03.73 

-1-31:16:03.8 

-1-6.0 

235 

0.34 

21 

Class 0/1 

9,10 

OMC-2 FIR 4 

05:35:26.97 

-05:09:54.5 

+11.4 

420 

30 

100 

IM proto-cluster 

11,12 

Cep E 

23:03:12.80 

-1-61:42:26.0 

-10.9 

730 

35 

100 

IM protostar 

13 

TIMASSS 

116293 

16:32:22.6 

-24:28:33 

-1-4.0 

120 

3 

22 

Class 0, HC 

14,15 


*PSC: pre-ste llar co r e; HC : h ot corino; WCCC: w ar m carbon-chain c hemis try ; IM: interm e diate- mass. 
References: 1197^) , ^|Evaris_et_alJ , ^Shirle^_et_aU (|200j 

^Marcelino et ah j 200^, jKrister^^^^Il ^ ^01 


^ T6t^_et_ay_^ll200^ ^ffirario_et_^U (199E), 
200^), ^‘^ Chen et~^ 1200E), 


^Grimier et alj 2009t) . ^ ^Furlan et alj (12014 ') . ^ ^Crimier et alj ll2010a G ^ iLoinard et alj (|2008| G ^ ^Correia et alj (120041 '). 


we detected transitions with upper level energies, Eup, be¬ 
low 150 K, and spontaneous emission coefficients, Aij, above 
10“® s“^ and 5 x 10“® s“^, respectively for the 2/3-mm and 
the 1-mm data. For HNCO, we detected transitions with 
Eu < 150 K and Aij > 10“® s“^. Tables iBl] and IB^ list 
all the NH 2 CHO and HNCO transitions fulfilling these cri¬ 
teria in the observed millimetre bands, as well as the 3 (t 
detections for each source. The sources where no NH 2 CHO 
lines were detected (see below) are not included in Table [Bl] 
For some sources with not many clear formamide detections 
(e.g. IRAS4A, Cep E), we included a few additional lines 
with peak intensities between 2a to 3c7, as indicated in the 
tables. We then htted the lines with a Gaussian function, 
and excluded from further analysis those falling well below 
or above the systemic velocity, and/or displaying too narrow 
or too broad line widths with respect to the typical values 
encountered for each source. 

Table [2] lists, for each source, the number of NH 2 CHO 
and HNCO lines detected and used in our analysis lSect l4.2Tl . 
While HNCO is easily detected in all the sources, NH 2 CHO 
remains undetected in five objects; L1544, TMC-1, Bl, 
L1527, and L1157mm. Moreover, in those sources where it is 
detected, the lines are typically weak (S/N~ 3 — 5). OMC- 
2 FIR 4 has the highest number of detected formamide lines, 
which are also the most intense. The results from the Gaus¬ 
sian fitting to the detected lines are presented in Tables IB3I 
- IB12I A sample of lines for all the ASAI sources are shown 
in Figs. ET]-El 


4.2 Derivation of physical properties 

Rotational diagram analysis 

In order to determine the excitation conditions -i.e. exci¬ 
tation temperature, column density and, eventually, abun¬ 
dance with respect to H 2 - of NH 2 GHO and HNGO for each 
source in a uniform way, we employed the GASSIS software 
to build rotational diagrams. This approach assumes (i) that 


Table 2. Number of NH2CHO and HNCO detected lines 


Source 

NH2CHO 
# Eu (K) 

HNCO 

# (k) 

L1544® 

0 

— 

2 

10-16 

TMCl 

0 

— 

3 

10-16 

Bl 

0 

— 

4 

10-30 

L1527 

0 

— 

4 

10-30 

L1157-mm 

0 

— 

4 

10-30 

IRAS 4A 

7 

15-70 

10 

10-130 

SVS13A 

13 

15-130 

19 

10-130 

OMC-2 FIR 4 

21 

10-130 

9 

10-100 

Cep E 

5 

10-22 

5 

10-85 

116293 

12 

10-160 

16 

10-95 


®Only 3-mm data available. 


the lines are optically thin, and (ii) Local Thermodynamic 
Equilibrium (LTE), meaning that a single Boltzmann tem¬ 
perature, known as rotational temperature, describes the rel¬ 
ative distribution of the population of all the energy levels 
for a given molecule. Under these assumptions, the upper- 
level column density 


Snku'^ 1 
hc^Aui ?7bf 


TnbdU 


( 1 ) 


and the rotational temperature. Trot, are related as follows: 

In ^ = ZniVtot - InQiE^t) - (2) 

Qu f^J-rot 

where k, u, h, and c are, respectively, Boltzmann’s constant, 
the frequency of the transition, Planck’s constant, and the 
speed of light; Qu is the degeneracy of the upper level, and 
Ntot is the total column density of the molecule. The second 
fraction in Eq. [T] is the inverse of the beam-hlling factor. 
We estimated it assuming sources with a gaussian intensity 
distribution: 
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el 

'nht Q2 Q2 (3) 

with 6a and 0b being, respectively, the source and telescope 
beam sizes. We adopted the source sizes indicated in Ta¬ 
ble O In those sources where a hot [T > 100 K) inner 
region is believed to exist, we considered two possible so¬ 
lutions: (i) the emission originates from a compact size rep¬ 
resenting this inner region or hot corino, which typically 
shows enhanced abundances of Complex Organic Molecules 
(COMs); and (ii) the emission homogeneously arises from 
the entire extended molecular envelope of the protostar. We 
determined the sizes of the compact hot corino regions ei¬ 
ther from published interferometric maps (SVS13A) or from 
the gas density structure, n(r), reported in the literature 
(116293, IRAS 4A, OMC-2, Cep E), as indicated in TablejS] 
In the latter case, we assumed a size equal to the diameter 
within which the dust temperature is above 100 K. 

Some sources, such as IRAS 4A, OMC-2, and Cep E, 
show extended velocity wings in a few of their lines. In order 
to separate their contribution to the line emission, we deter¬ 
mined their line flux, J T^hdV, by fitting a Gaussian func¬ 
tion to the affected lines after masking their high-velocity 
wings. In sources with two to four well-aligned data points 
in the rotational diagrams, we took into account the rela¬ 
tively large error bars by fitting two additional “extreme” 
lines passing through the tips of the error bars of the lowest 
and largest energy points. An example is shown for B1 in 
Fig.[T] where the two extreme solutions are depicted in blue, 
while the best solution is marked in red. The remaining ro¬ 
tational diagrams and the best fit to their data points using 
Eq. [2] are shown in Figs. IC4I and IC5I where the error bars 
take into account calibration errors as well as the rms value 
around each line. 

We co mpared our rotat i onal d iagram results with those 
reported in iMarcelino et ^ ll2009l ~) for the four sources com¬ 
mon to both studies: L1544, TMC-1, Bl, and L1527. The 
column densities of HNCO are in perfect agreement within 
the uncertainties, while the rotational temperatures agree 
within 1 K. 

For homogeneity with the methodology used for 
NH 2 CHO, we estimated the properties of HNCO in the LTE 
approximation. In addition, by adopting the same source 
sizes for HNCO and NH 2 CHO, we assumed that the emis¬ 
sion from both molecules originates in the same region(s). 
The similar average line widths between the two species sug¬ 
gest this is a reasonable assumption. Tableland Figs. IC4I 
and IC5I present the results of the rotational diagram anal¬ 
ysis. For most of the sources, a single component fits well 
both the NH 2 CHO and HNCO points and therefore LTE 
seems to reproduce well the observations. This can also be 
seen in Figs. lCll to lC3l where the observed spectra (in black) 
and the best fit models (in red) match fairly well. However, 
for SVS13A, Cep E, and OMC-2, the compact solutions cor¬ 
respond to HNCO lines that are moderately optically thick 
(r ~ 1 —10). The most extreme case is Cep E, for which also 
the NH 2 CHO lines are optically thick. This is in contradic¬ 
tion with the underlying assumption of optically thin lines 
in the rotational diagram method. We find, however, that 
this caveat can be easily overcome by adopting a slightly 
larger source size, of 3”, 2", and 2", respectively for SVS13A, 
Cep E, and OMC-2. Doing this, the resulting column den- 



10 20 30 


Eu/k (K) 

Figure 1. HNCO rotational diagram of Bl. Data points are de¬ 
picted in black. The red lines correspond to the best fit to the 
data points. The extreme solutions taking into account the error 
bars are displayed in dashed blue. 


sities are reduced by a factor 2 (OMC-2) to 15 (Cep E), r 
becomes much smaller than 1, and the lines can be well fit¬ 
ted by the solutions. Consequently, the uncertainties in the 
compact-solution column densities in these three sources are 
larger than reported in Table |3l but they are taken into ac¬ 
count in the discussion (Sect.[S] see Figs. [5] and O. 

There are two objects where a single component does 
not appear to explain the emission of all the lines: IRAS 4A 
and 116293, two well-known hot corino sources. Indeed, their 
rotational diagrams suggest either the contribution of two 
components, or non-LTE effects, coming into play. Con¬ 
sidering the former. Table |3] presents the results of a two- 
component solution to the rotational diagrams of these two 
objects, where Cl is assumed to represent the cold extended 
envelope of the protostar, and C2 the small inner hot corino. 
While this 2-component solution reproduces well the obser¬ 
vations, non-LTE effects cannot be ruled out. 

As for the five objects where formamide was not de¬ 
tected, we determined a 3cr upper limit to its column density 
under the assumption of LTE and adopting the correspond¬ 
ing value of Trot derived for HNCO. To this end, we used 
the spectral data around the NH 2 CHO 4o,4 - 3o,3 transi¬ 
tion at 84.542 GHz, expected to be the most intense at the 
cold temperatures implied by the HNCO results. The upper 
limits thus derived are shown in Table |3] 

Once the column densities of HNCO and NH 2 CHO were 
obtained, we derived their respective abundances with re¬ 
spect to molecular hydrogen (H 2 ) using the H 2 column den¬ 
sities, A^h 2 , listed in Table [3] which correspond to the indi¬ 
cated source sizes. The uncertainty on Ah 2 is included for 
those sources where this was provided in the corresponding 
bibliographic reference. The resulting abundances span more 
than two orders of magnitude and are shown in Table |3l to¬ 
gether with their ratio, R = X(HNCO)/X(NH 2 CHO). 
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Table 3 . Results from the rotational diagram analysis of NH 2 CHO and HNCO: Adopted size and H 2 column densities (A^h2)5 derived 
rotational temperatures, Trot, derived HNCO and HN 2 CHO column densities (A^hnco? ^nh2CHo)? resulting abundances with respect 
to H 2 ( A’hnco, A’nh2CHo), ^rid ratio of HNCO to NH 2 CHO column densities (R). 


Source 

Size® 

(") 

AJO 

-^'^H2 

(10^^cm~^) 

Erot(HNCO) 

(K) 

AhnCO 

(lO^^cm”®) 

AhnCO 

(10-11) 

rnc 

-^rot(NH2CHO) 

(K) 

^'^NH2CHO 

(10^^cm“^) 

AnH2CHO 

(10-11) 

R 

1-component fit 

L1544 

BF 

9.4± 1.6i 

7± 3 

5±3 

5±3 

7 

< 0.036 

< 0.046 

> 130 

TMCl'i 

BF 

1.0 ±0.1^ 

4± 1 

8±5 

80± 50 

4 

< 0.47 

< 5.2 

> 17 

B1 

BF 

7.9 ±0.3® 

10 ±2 

8.4 ± 1.6 

11 ± 2 

10 

< 0.087 

< 0.11 

> 97 

L1527 

BF 

d.D 

7.5 ± 1.4 

2.5 ± 1.5 

6±4 

7.5 

< 0.062 

< 0.15 

> 40 

L1157-mm 

30 

120® 

8± 1 

4± 1 

0.35 ± 0.03 

8 

< 1 

< 0.008 

> 40 

SVS13A (ext)® 

20 

10® 

58 ±6 

11 ± 2 

11 ± 2 

64 ±6 

3.0 ±0.4 

3.0 ±0.4 

4± 1 

(com)® 

1 

1000'^ 

36 ± 3 

1500 ± 300 

15 ±3 

40 ±4 

320 ± 60 

3.2 ±0.6 

5± 1 

OMC-2 (ext) 

25 

19® 

25 ± 3 

16 ±3 

1.9 ± 0.4 

58 ±4 

3.1 ±0.2 

0.36 ±0.02 

5± 1 

(com) 

2 

4.6® 

19 ± 1 

900 ± 100 

910 ± 80 

32 ±2 

110 ± 10 

110 ± 10 

8± 1 

Cep E (ext) 

40 

4.8^ 

30 ±5 

6.2 ±0.3 

13 ± 1 

9±2 

0.2 ±0.1 

0.4 ±0.2 

30± 13 

(com) 

0.5 

230® 

17 ± 1 

6000 ± 1000 

130 ± 15 

6± 1 

500 ± 300 

11±5 

12 ±6 


2-component fit 


IRAS 4A (Cl) 

30 

2 . 91 ° 

11 ±3 

10 ± 1 

IRAS 4A (C2) 

0.5 

2501® 

43 ±8 

2000 ± 1000 

116293 (Cl) 

30 

2 . 9 II 

14 ±5 

20 ±2 

116293 (C2) 

1.2 

53 II 

47 ±4 

4400 ± 700 


34 ±2 

19± 15 

0.6 ±0.5 

1.9 ±0.2 

18 ±2 

80 ±40 

30 ±5 

500 ± 100 

20 ±5 

4±2 

69 ± 7 

5± 1 

1.7 ±0.6 

6±2 

12 ±4 

830 ±130 

83 ±33 

590 ± 190 

110 ±40 

8±3 


^BF: beam-filling assu med. For the other source s, the size has been adopted as follows: L115 7mm and IRAS 4A (extended) from 
Ijorgensen et al] j^OO^; SVS1 3 A (ex tended) from iLefloch et alJ lll99^: SVS13A (compact ) from [Looney et al.lll20 QOlh OM C-2 FIR 4 
(extended) fror n iFurlan et ^ |j[20]^; OMC-2 FIR 4 ( compact) from [Crimie r et alJ ([200^; Cep E from lchjnTrt alJ (l20oi]) : IRAS 4A 
(compact) fromjMare^^an~ gOod): IRAS 16293 from [jaberet al] l|2014|). 

^References: ^CmpstT^S] 120^), ^Maezawa et_aJj 1^99^. y3ame]_eCalJ ll20lj), ^Parise et al.| (|2014[). ^Jfirpjensen et ahj l|2002|) 
^Lefloch et alJ jl998l ). iLoonev et alJ ll200nh . ^Crimie^r et al.l ll2009h . ^Crimier et al.l l^2010a^ . ^ ^Maret et alJ ll2002l ). ^ ^Crimier et all 

2010hl) . 


"^For the non-detections of NH 2 CHO, we have computed a Scr upper limit to its column density adopting the same Trot derived for 
HNCO (see text). 

“^Data for NH 2 CHO upper limit derived from 3-mm data by N. Marcelino. 

^^^(HNCO) is probably a lower limit due to contamination from the OFF position. 


4-2.2 Radiative transfer analysis taking into account the 
source structure 

The source structures of 11 6293, IRAS 4A, and OMC-2, are 
report ed in the literature (iMaret et al.ll20(d . ICrimier et al.l 
l2010a ). Therefore, for these objects, a more sophisticated ra¬ 
diative transfer analysis is possible that takes into account 
the temperature and gas density as a function of distance 
to the central protost ar. Cep E also has a known structure 
llCrimier et al.l BoiOal ). but having only 5 line detections, 
both in HNCO and NH 2 CHO, we do not consider it here. 
We have analysed the 116293, IRAS 4A, and OMC-2 lines 
by means of the code GRAPES (GRenoble Analysis of Pro- 
tostellar Envelope Spe ctra) , whose details are described in 
ICeccarelli et al.l ll2003l ) and (jaber et al.l (l2014l ). 

Briefly, GRAPES computes the Spectral Line Energy 
Distribution (SLED) of a free-infalling spherical envelope 
with given gas and dust density and temperature profiles, 
and for a given mass of the central object. The dust to gas 
ratio is assumed to be the standard one, 0.01 in mass, and 
the grains have an average diameter of 0.1 ^m. The species 
abundance is assumed to follow a step-function, with a jump 
at the dust temperature Tjump, which simul ates the thermal 
desor ption of species from icy mantles (e.g. ICeccarelli et al.l 
I 2 OOOI ). The abundance Xi in the warm (T ^ Ijump) enve¬ 
lope is constant. In the outer envelope, we assumed that the 


abundance follows a power law as a function o f the radius, 
AotL with an index equal to 0, -1 and -2, as in ijaber et al.l 
ll2014l) . Ai and Ao are considered parameters of the model. 
Since, to our knowledge, the binding energy of NH 2 CHO is 
not available in the literature, we treat Tjump as a parameter 
too. However, if the molecules are trapped in water ice, the 
binding energy of H 2 O will largely determine the dust tem¬ 
perature at which NH 2 GHO is injected into the gas phase. 

The radiative transfer is solved with the escape prob¬ 
ability formalism and the escape probability is computed 
integrating each line opacity over the dvr solid angle. We ran 
models assuming LTE populations for formamide and, for 
comparison with Sect. 14.2711 HNGO, and models taking into 
account non-LTE effects for HNG O. In the latt er case, we 
used the collisional coefficie nts by &een| Il98^, retrieved 
from the LAMDA database (ISchoier et al.lf^05l ). 


For each molecule and source, we ran a large grid of 
models varying the four parameters above: Ai, Ao, Ijump, 
and a. In total, we ran about 20,000 models per source. The 
computed SLED of each model was then compared with the 
observed SLED to find the solution with the best fit. The 
results of this analysis are reported in Table [d] where we 
give the best fit values and the range of Ai, Ao , Tjump with 
^ 1. We note that the there is no appreciable difference 
in the best when using a different value of a, so we took 
the simplest solution: a = 0. In this respect, the situation is 
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Table 4. Results of grapes analysis for NH 2 CHO and HNCO 
considering the source structure of IRAS 4A, 116293 and OMC2* 



IRAS 4A 

116293 

OMC-2 

HNCO LTE 

Ao (10-11) 

3± 1 

0.1 ± 0.1 

5.5 ± 1.5 

A; (10-11) 

20 ± 10 

90± 10 

< 170 

Tjump (K) 

100 

40 

80 

Tjump range (K) 

60 - 120 

30 - 50 

Ss 30 


1.2 

2.0 

1.0 

HNCO non-LTE 

Ao (10-11) 

3± 1 

0.5 ± 0.4 

4± 1 

A; (10-11) 

30 ±20 

600 ± 300 

< 20 

Tjump (K) 

100 

90 

80 

Tjump range (K) 

^ 50 

^ 60 

Ss 30 

X^ 

1.0 

1.5 

0.7 

NH 2 CHO 

Ao (10-11) 

2 ± 1 

0.3 ± 0.2 

0.3 ± 0.3 

A; (10-11) 

50± 10 

60 ±20 

200 ± 50 

Tjump (K) 

100 

90 

80 

Tjump range (K) 

> 100 

Ss 50 

60 - 100 

X^ 

2.0 

0.7 

1.3 

R = 

A(HNCO)/A(NH2CHO) 


Ro (LTE) 

1.5 ±0.9 

< 1.7 

18± 18 

Ri (LTE) 

0.4 ±0.2 

1.5 ± 0.5 

< 0.85 

Ro (non-LTE) 

1.5 ±0.9 

1.7 ± 1.7 

13 ± 13 

Ri (non-LTE) 

0.6 ±0.4 

10 ±6 

< 0.1 


‘Abundances with respect to H 2 are times 10 Xo and A; 
are the outer and inner abundances, respectively. 

similar to what Ijaber et al.l (l2014n found in their study of 
IRAS16293. 

A comparison between the results obtained for HNCO 
with the LTE and non-LTE level populations shows that the 
LTE approximation is quite good in the case of IRAS 4A and 
OMC-2, but not for 116293. The reason for that is probably 
a lower density envelope of 116293 compared to the other 
two sources. Therefore, the LTE results are likely reliable 
also for the formamide in IRAS 4A and OMC-2, while in 
116293 these have to be taken with some more caution. 

A second result of the GRAPES analysis is that both 
HNCO and formamide have a jump in their abundances at 
roughly the same dust temperature, 80-100 K. This is an 
important result reflecting the two molecules have similar 
behaviours with changes in temperature. It suggests they 
trace the same regions within the analysed protostars. 

In order to evaluate whether the rotational diagram 
(hereafter RD) and GRAPES analyses are in agreement, we 
compare their respective abundance values, which roughly 
agree within an order of magnitude, in the Appendix|^ We 
note here that, while the GRAPES analysis is likely more 
accurate, we are not able to apply it to the other sources of 
this study due to the lack of known source structure and/or 
lack of a sufficient amount of molecular lines. The absence of 
interferometric imaging of the HNCO and NH 2 CHO emis¬ 
sion also hinders the study of the inner structure of the pro- 
tostellar emission. Therefore, we base the discussion below 
largely on the RD results, with a note of caution that those 


values may not strictly represent the physical properties of 
the sources. 

5 DISCUSSION 

5.1 Formation routes of NH 2 CHO 

The formation mechanism(s) of interstellar formamide, as 
that of other COMs, is still far from being established. Sev¬ 
eral routes have been proposed so far which include both 
gas- phase and grain- s urface processes. Concerning the for¬ 
mer, lOuan fc HerbstI |2003) suggested NH 2 CHO forms via 
the radiative association reaction 

NH+ -b H2CO ^ NH4CH2O+ -b hi/ ( 4 ) 

followed by dissociative recombination. iHalfen et al.l (1201 ifl 
proposed the following ion-molecule reaction and subsequent 
electron recombination: 

NH+ -b H 2 CO ^ NH 3 CHO+ -b H 2 (5) 

NH 3 CHO+ -b e" ^ NH 2 CHO -b H ( 6 ) 

These reactions all have unknown rates. Thus, further ex¬ 
perimental work will be needed in order to evaluate their 
effectiveness in producing formamide. 

Neutral neutral reactions have also been discussed as 
possible gas-phase r outes leading to NH 2 CHO. In particular, 
iGarrod et all ll2008ll proposed the radical-neutral reaction 

H 2 CO -b NH 2 ^ NH 2 CHO -b H (7) 

However, as recently mentioned bv iRedondo et ^ ll2014h . 
it presents a net activation barrier of > 1000 K that makes 
it inviable in interstellar conditions. Other neutral-neutral 
reactions evaluated by these authors also revealed to have 
large activation barriers, thus ruling them out as dominant 
or efficient mechanisms to produce NH 2 GHO. 

Eorma mide may also be f ormed on the icy mantles of 
dust grains. I Jones et al.l ll201ll ') conducted some experimen¬ 
tal work in which they irradiate a mixture of ammonia 
(NH 3 ) and carbon monoxide (GO) ices with high-energy 
(keV) electrons, resulting in NH 2 GHO as one of the final 
products. The authors discuss several possible reactions and 
conclude that the most plausible route towards formamide 
begins with the cleavage of the nitrogen-hydrogen bond of 
ammonia, forming the NH 2 radical and atomic H. The lat¬ 
ter, containing excess kinetic energy, can then add to CO, 
overcoming the entrance barrier, to produce the formyl rad¬ 
ical (HCO). Einally, HCO can combine with NH 2 to yield 
NH 2 CHO. 

A different gra in-mantle mechanism was proposed by 
iGarrod et ahl (l2008ll . who considered hydrogenation (i.e. ad¬ 
dition of H atoms) of OCN in their chemical models. How¬ 
ever, this route resulted in an overabundance of NH 2 GHO 
and an underabundance HNGO, since the latter was effi- 
ciently hydrogenated to yield formamide, the final product. 
iRaunier et ahl ||2004|') performed experimental Vacuum Ultra 
Violet (VUV) irradiation of solid HNCO at 10 K, which led 
to NH 2 CHO among the final products. They proposed that 
photodissociation of HNCO yields free H atoms that subse¬ 
quently hydrogenate other HNCO molecules in the solid to 
finally give NH 2 CHO. The limitation of this e xperiment is 
that it was carried out with pure solid HNCO. Ijones et ahl 
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1 I 2011 I) mentioned that, in the presence of NH 3 , quite abun¬ 
dant in grain mantles, HNCO will preferentially react with 
it, resulting in NH^-f OCN“. Despite these caveats, hydro¬ 
genation of HNCO on grain mantles was recently found to 
be a most likely solution i n the case of the ou tflow shock 
regions L1157-B1 and B2 (iMendoza et aLl|2014ll . More ex¬ 
periments and calculations are needed in order to assess the 
efficiency of this formation route. 


5.2 Correlation between HNCO and NH 2 CHO 


From the previous section, it is clear that, until more gas- 
phase and surface reaction rates involving the mentioned 
species are measured, it will be difficult to establish the exact 
synthesis mechanisms of formamide in space. 

In this section, we assess, from an observational point of 
view, whether hydrogenation of HNCO leading to NH 2 CHO 
on the icy mantles of dust grains could be a dominant for¬ 
mation route. To this aim, we plot in Fig. [2] the NH 2 CHO 
versus HNCO abundances of all ou r sources (Table [3l), a s 
well as the shock regions analysed by[Mendoza^_et_^ (uO^), 
and the high-mass SFRs re ported in iBisschop et ail (l2007l l 
and iNummelin et al.l (l2000l i. for comparison. The latter were 
obtained by the cited authors via the RD method assuming 
the emission comes from the inner hot core regions. Thus, 
for homogeneity, we split the plot into two panels, the up¬ 
per one showing only the compact/inner solutions of the 
RD analysis, classified by masses. The best power law ht 
to these points is marked with a dashed line, and is given 
by the equation X(NH 2 CHO) = 0.04X(HNCO)°■®^ with 
a Pearson coefficient of 0.96, indicating a tight correlation. 
The fact that this correlation is almost linear and holds for 
more than three orders of magnitude in abundance suggests 
that HNCO and NH 2 CHO are chemically related. This re¬ 
sult confi rms, on a more s tatist ical basis, what was recently 
found by IMendoza et ^ (120141 '). 

However, this correlation does not hold for the ob¬ 
jects without formamide detections, which are plotted in the 
lower panel of Fig. [5] together with the extended envelope 
solutions of the RD analysis. Here, it is clearly seen that 
all the upper limits lie well below the best fit line, indicat¬ 
ing a significantly larger amount of gas-phase HNCO rela¬ 
tive to NH 2 CHO in comparison to the other sources. These 
objects are the coldest in our sample, representing either 
pre-stellar cores or protostars with no detectable hot corino 
within them. The rotational temperatures inferred from the 
HNCO RD analysis are also among the lowest in our sample. 
In this same plot, the points representing formamide detec¬ 
tions (extended envelope component) also show a tendency 
towards lower relative values of X(NH 2 CHO), although not 
as pronounced. 

Thus, it appears that regions with colder temperatures 
are more dehcient in NH 2 CHO than protostars with hot in¬ 
ner regions, indicating that higher temperatures are needed 
for NH 2 CHO to become relatively abundant in the gas 
phase. This might be explained by (i) NH 2 CHO forming 
in the gas phase at temperatures above 100 K, and/or (ii) 
it forming predominantly on the icy mantles of dust grains 
at low temperatures, and subsequently sublimating into the 
gas-phase when the temperature in the inn er regions rises 
suffici ently. As for the former possibility, IMendoza et al.l 
(l2014h quantitatively argued that reaction [7] does not suf- 



X(HNCO) 

Figure 2. Plot of NH 2 CHO versus HNCO abundances with 
respect to H 2 . Top-. Data points included in the power-law fit 
(dashed line; see text). Red squares and green diamonds denote 
the compact or inner RD solutions of low- and intermediate-mass 
sources in this study, respectively. Magenta triangles and black 
stars correspond, re spectively, to outflow shock regions (from 
Mendozaet_^ I 2 OI 4 II and hi^h-m ass sources ('from iBisschop et alj 
20071 and iNummelin et al.ll2000h . Bottom: Data points not in¬ 
cluded in the power-law fit (see text). Blue open squares represent 
the extended or outer RD solutions, while black open and filled 
circles denote the GRAPES LTE values for the outer and inner 
components, respectively. 


fice to explain the amount of gas-phase formamide in the 
shock regions of L1157 protostellar outflow. In addition, 
the high activation barrier the reaction needs to overcome 
makes this an unviable route. Other purely gas-phase for¬ 
mation routes still need more investigation in terms of reac- 
tion rates and a ctivat ion barriers, as discussed in Sect. o 
IMendoza et al.l ll2014ll favoured a grain formation mecha¬ 
nism followed by mantle-grain evaporation/spnttering on 
the basis of the comparable abundance enhancements of 
HNCO, NH 2 CHO, and CH 3 OH in the gas-phase between 
the two protostellar shocks studied by the authors. There- 
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fore, grain formation of NH 2 CHO appears to be the most 
likely possibility. 

On the other hand, while grain formation of HNCO 
is likely to occur in the cold phases of star formation 
llHaseeawa fc HerbstlllQQ^l . gas-phase reactions leading to 
HNCO at such col d temperatures can also take place 
efficiently (see e.g. iMarcelino et al.l l2009l and references 
therein), overcoming strong depletion. This would explain 
its relatively high gas-phase abundance already in the very 
early -and cold- phases of star formation, and also the high 
values of HNCO to NH 2 CHO abundance ratios we find in 
the coldest sources of our sample. 

In Fig. [3 we plot the HNCO abundance, the NH 2 CHO 
abundance, and their ratio, R, as a function of bolomet- 
ric luminosity for those sources with a reported luminosity 
estimate (see Table [1]). For the objects in our study with 
formamide detection, we only plot the points correspond¬ 
ing to the inner or compact component (red circles), since 
these regions are expected to be the dominant contributors 
to the overall luminosity. The HNCO and NH 2 CHO abun¬ 
dance panels both show the high-mass sources lying on top 
of the plot, while the points representing our sample sources 
are more scattered, with the coldest objects (in blue) show¬ 
ing the lowest abundances. This trend is much more pro¬ 
nounced in the case of NH 2 CHO, for which hot corino re¬ 
gions (red points) display higher NH 2 CHO abundances than 
the colder objects by more than an order of magnitude. More 
interesting is the plot of R, which illustrates how this quan¬ 
tity remains roughly constant along 6 orders of magnitude 
in luminosity for the NH 2 CHO-emitting sources, with values 
ranging from 3 to 10 approximately. This reflects the almost- 
linearity of the correlation between the abundance of the two 
species. On the other hand, this value rises considerably for 
the lower luminosity sources, re-enforcing our interpretation 
that formamide mostly forms on grains at cold tempera¬ 
tures, while HNCO may form both on grains and in the gas. 

The strikingly tight and almost linear correlation be¬ 
tween the abundance of the two molecules once NH 2 CHO 
becomes detectable suggests one of the two following possi¬ 
bilities: (i) HNCO and NH 2 CHO are both formed from the 
same parent species on dust grain mantles, or (ii) one forms 
from the other. Among the grain formation routes that have 
been proposed so far, hydrogenation of HNCO leading to 
NH 2 CHO is the only mechanism that would explain our ob¬ 
servational results. While this route is found to have some 
caveats (see Sect. 15.iT) . it is also true that more experimental 
work is needed to better assess its efficiency. 

If the abundance of gaseous NH 2 CHO truly depends on 
temperature, we should find a difference in R between the 
hot corino and the cold envelope regions of IRAS 4A and 
116293. Looking at Table O this is indeed the case. As for 
OMC-2, Cep E, and SVS13A, only one component was nec¬ 
essary to describe their rotational temperatures and column 
densities. Therefore, we cannot compare the extended and 
compact values as in the case of a 2-component solution. We 
can nevertheless guess that, excluding the case of Cep E, for 
which only low-energy formamide lines were detected, the 
compact solution is likely the best, given the low values of 
R and the relatively high rotational temperatures derived. 
This would imply most of the emission arises in the inner 
hot corino regions. For Cep E, more molecular observations 
at higher frequencies are needed to confirm this. 



Figure 3. Abundance of HNCO (top), NH 2 CHO (middle) and 
their ratio (bottom) against bolometric luminosity. Symbols are 
as in Fig. [2] 
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Figures [2] and [3] also include the results from the LTE 
GRAPES analysis. 116293 (labelled in the plots) is included 
for completeness despite the fact that the GRAPES analysis 
suggests non-LTE effects should be taken into account for 
this object. While these points introduce more scatter in 
the plots, it can be clearly seen that the inner components 
of the sources analysed with GRAPES have a lower HNGO 
abundance relative to NH 2 CHO, compared to what is found 
via the RD analysis. This yields lower R values, indicating a 
considerable amount of formamide with respect to HNGO in 
these regions and suggesting, as mentioned in Sect. m that 
the 2-component approximation in the RD analysis is over 
simplistic: while we assumed that only the higher-energy 
formamide lines arose from the compact inner region, it is 
likely that a significant amount of emission from the low- 
energy lines also originates here and not exclusively in the 
outer envelope. 

The trend showing higher R in the outer envelope than 
in the inner regions holds for both IRAS 4A and OMC- 
2, which further supports the fact that NH 2 GHO requires 
higher temperatures than HNGO to be detectable in the 
gas phase. This kind of analysis, taking into account the 
source structure, is needed in a larger sample of objects in 
order to draw conclusions about both the chemistry and 
the validity of our RD analysis on a more statistical basis. 
Interferometric mapping would also greatly help disentan¬ 
gling source multiplicity and verifying whether the emission 
of HNGO and NH 2 CHO trace the same regions, as has been 
assumed in this work. 


6 CONCLUSIONS 

As part of the IRAM Large programme ASAI, we searched 
for millimetre spectral lines from formamide (NH 2 CHO), a 
presumably crucial precursor of pre-biotic material, and iso- 
cyanic acid (HNGO), in ten low- and intermediate-mass star 
forming regions with different properties. The dataset, ob¬ 
tained with the IRAM 30-m telescope, consists mainly of 
unbiased broadband spectral surveys at 1, 2, and 3-mm. 
Our aim was to investigate the chemical connection be¬ 
tween these two molecular species and gain some observa¬ 
tional insights into the formation mechanisms of formamide 
in interstellar conditions. The present work represents the 
first systematic study within ASAI and statistically com¬ 
pletes the low-mass end of similar studies performed towards 
high-mass star-forming regions. Our main findings are sum¬ 
marised as follows. 

1. The high sensitivity and large frequency range of the 
spectral surveys allowed us to evaluate the detectability of 
numerous NH 2 GHO and HNGO transitions. We detect for¬ 
mamide in five out of the ten objects under study (IRAS 4A, 
IRAS 16293, SVS13A, Gep E, and OMG-2), and HNGO in 
all of them. Since formamide had already been detected in 
IRAS 16293 -also investigated here for completeness-, this 
study raises the number of known low- and intermediate- 
mass formamide-emitting protostars to five, thus signifi¬ 
cantly improving the statistics. 

2. We derived HNGO and NH 2 GHO column densities via 

the rotational diagram method for all the sources. As a re¬ 
sult, we found NH 2 GHO abundances with respect to H 2 in 
the range - 10“®, and HNGO abundances between 


and 10“®. For those objects without formamide de¬ 
tection, we provided an upper limit to its column density 
and abundance. 

3. For three targets (IRAS 4A, IRAS 16293, and OMG- 
2 ), the source density and temperature structures are known 
and published, and we were thus able to take them into ac¬ 
count through a more sophisticated analysis using the code 
GRAPES. This method fits an abundance profile that con¬ 
sists of a step function, with the separation between the 
two values roughly corresponding to the hot corino size. A 
comparison between the two radiative transfer analyses em¬ 
ployed reveals overall agreement within an order of magni¬ 
tude. The GRAPES analysis also indicates that one of the 
studied objects, IRAS 16293, requires a non-LTE radiative 
transfer analysis, which at the moment is not possible due 
to the lack of collisional coefficients for NH 2 GHO. LTE ap¬ 
pears to describe correctly the other two sources analysed 
with GRAPES, and is assumed to be a good approximation 
for all the other sources in our sample. 

4. For the sources where formamide was detected, i.e. 
hot corino sources, we found an almost linear correlation 
between HNGO and NH 2 GHO abundances that holds for 
several orders of magnitude. This suggests that the two 
molecules may be chemically associated. On the other hand, 
those sources with no formamide detection do not follow this 
correlation, but instead show much larger amounts of HNGO 
relative to NH 2 GHO. These objects are the coldest in this 
study, and unlike the rest of our sample, they contain no 
known hot corinos. 

5. Our findings and the NH 2 GHO formation routes pro¬ 
posed so far in the literature suggest that, unlike HNGO, 
NH 2 GHO does not form efficiently in the gas phase at cold 
temperatures and may be formed on the mantles of dust 
grains, where it remains frozen at cold temperatures. As 
soon as the temperature rises sufficiently to sublimate the 
icy grain mantles, formamide is incorporated into the gas 
and becomes detectable. The tight and almost linear cor¬ 
relation with HNGO suggests a possible formation route of 
NH 2 GHO via hydrogenation of HNGO, although other pos¬ 
sibilities should not be ruled out. In particular, two poten¬ 
tially viable gas-phase pathways leading to formamide in¬ 
volve formaldehyde (H 2 GO). It is therefore worth exploring 
the connection between H 2 GO and NH 2 GHO, which will be 
the subject of a forthcoming paper. 

6 . In order to evaluate the validity of our conclusions, sev¬ 
eral aspects need to be explored more thoroughly. From an 
observational point of view, interferometric imaging is nec¬ 
essary to assess the relative spatial distribution of HNGO 
and NH 2 GHO, and retrieve more accurate abundance ra¬ 
tios, in particular in the hot corino sources. In addition, 
more detailed and sophisticated radiative transfer analysis 
requires, on the one hand, knowledge of the source density 
and temperature profiles, and on the other hand, collisional 
coefficient calculations for NH 2 GHO, currently unavailable. 
Finally, more chemical experiments are needed to estimate 
the efficiency of the hydrogenation processes leading from 
isocyanic acid to formamide on interstellar dust grains, as 
well as the viability of purely gas-phase reactions. 
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Table Al. Comparison between GRAPES and RD analyses* 



IRAS 4A 

OMC-2 

Inner size (") 

1.5 

3.1 

RD-to-GRAPES ratio (LTE) 

a:o(hnco) 

11 ±4 

1.5 ±0.5 

A:i(HNCO) 

4±3 

> 2.5 

Ao(NH2 CHO) 

1±1 

5±5 

Ai(NH2CHO) 

0.4 ±0.4 

0.3 ±0.3 

Ro 

12 ±7 

0.3 ±0.3 

Ri 

11±8 

> 10 


RD-to-GRAPES ratio (non-LTE) 


a:o(hnco) 

11±4 

2.1 ±0.7 

A:i(HNCO) 

3±2 

> 22 

Ao(NH2 CHO) 

1±1 

5±5 

Ai(NH2CHO) 

0.4 ±0.4 

0.3 ±0.3 

Ro 

12 ±7 

0.4 ±0.4 

Ri 

7±6 

> 85 


‘Abundances with respect to H 2 are times 10 


APPENDIX A: COMPARISON BETWEEN 
GRAPES AND ROTATIONAL DIAGRAM 
ANALYSES 

This section aims to compare the agreement between the 
rotational diagram and GRAPES methods. As described in 
Sect l4.2.'2l the line emission in 116293 does not appear to be 
well described by LTE, and a more realistic radiative trans¬ 
fer treatment will need to wait until collisional coefficients 
are available for NH 2 CHO. Therefore, we do not consider 
it here, while it is worth noticing that a rotational diagram 
analysis is likely too simplistic to analyse the HNCO and 
NH 2 CHO lines in this source. 

In OMC-2, the GRAPES analysis tells us that the tem¬ 
perature that separates the inner and outer components is 
80 K, both in the LTE and non-LTE approximations. Thus, 
for consistency in the comparison, we re-computed the inner 
and outer abundances resulting from the RD analysis using 
the same inner sizes as in GRAPES, instead of those cor¬ 
responding to a temperature of the 100 K (see Sect. l4.2Tni . 
We note that, while the RD analysis allowed for a separation 
of two components (inner and outer) for IRAS 4A, a single 
component was sufficient for OMC-2. It should be kept in 
mind, therefore, that for the latter the comparison is not 
equivalent, since we are not comparing a two-component so¬ 
lution with another two-component solution as in the case 
of the other two protostars. 

The results of the comparison are listed in Table [AT] and 
illustrated in Fig. iah where we present the comparison us¬ 
ing both the LTE and non-LTE results from GRAPES. It is 
evident that LTE and non-LTE yield practically the same re¬ 
sults for these two sources. It can also be seen that the errors 
are quite high in some cases, up to 100%, which are caused 
by the large uncertainties resulting from the GRAPES anal¬ 
ysis. Taking these into account, we find the following be¬ 
haviours: 

• HNCO abundance: Generally, both methods agree 
within an order of magnitude, but there is a tendency to- 


RD—to—GRAPES ratio 


100 

■q 

y 10 


X 1 

0.1 


100 

o' 

S '0 

CM 

i 1 
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0.1 


100 


10 

cn 

1 


0.1 











1 . 0 '" 






Figure Al. Ratio of RD-to-GRAPES abundances. Top-. HNGO 
abundance. Middle: NH 2 GHO abundance. Bottom: HNCO to 
NH 2 CHO abundance ratio. Filled and open circles represent, re¬ 
spectively, the LTE and non-LTE HNCO solution in the GRAPES 
analysis. The horizontal dashed lines mark equality between RD 
and GRAPES values. 


wards higher values in the RD analysis, by a factor of a 
few. 

• NH 2 CHO abundance: Again, we find agreement within 
a factor of a few. The compact solution is systematically 
lower in the RD treatment. This suggests that a non- 
negligible amount of emission from low-energy molecular 
lines actually comes from the inner region, and not exclu¬ 
sively from the extended envelope, as assumed in the linear 
fitting of the RD. Such a finding reflects the necessity of 
analysis like that performed with GRAPES if we want to 
properly disentangle the inner and outer components in hot 
corino or hot core sources. 

• HNGO to NH 2 GHO abundance ratio, R: In this case, 
the two analysis methods agree within a factor of a few. 
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APPENDIX B: TABLES 



Formamide 


Table Bl. NH 2 CHO transitions searched for in this study and 3cr detections^ 


Transition 

V 

(MHz) 

Eu 

(K) 

Aui 

(10-® s-i) 

dh 

n 

OMC-2 

CepE 

SVS13A 

IRAS4A 

116293 

Blends 

4i.4 - 3i,3 

81693.446 

12.8 

3.46 

30 

N 

Y* 

N 

N 

N 


4o,4 - 3o,3 

84542.330 

10.2 

4.09 

29 

Y 

Y 

N 

N 

Y 


42.3 - 32,2 

84807.795 

22.1 

3.09 

29 

N 

N 

N 

N 

N 


43,2 - 33.1 

84888.994 

37.0 

1.81 

29 

N 

N 

N 

N 

N 


43,1 - 33.0 

84890.987 

37.0 

1.81 

29 

N 

N 

N 

N 

N 


42,2 - 32,1 

85093.272 

22.1 

3.13 

29 

N 

Y 

N 

N 

N 


4l.3 - 3 i,2 

87848.873 

13.5 

4.30 

28 

Y 

Y* 

N 

N 

Y 


5i.5 - 4i_4 

102064.267 

17.7 

7.06 

24 

Y 

Y 

N 

B 

Y 

H 2 C 0 H+? 

5o,5 - 4o,4 

105464.219 

15.2 

8.11 

23 

Y 

Y 

Y 

Y 

Y 


52,4 - 42,3 

105972.599 

27.2 

6.92 

23 

Y 

N 

N 

N 

N 


54,2 - 44,1 

106107.870 

63.0 

2.98 

23 

N 

N 

B 

B 

N 

NH 2 CH 0 

54,1 - 44,0 

106107.895 

63.0 

2.98 

23 

N 

N 

B 

B 

N 

NH 2 CH 0 

53,3 - 43_2 

106134.427 

42.1 

5.29 

23 

N 

N 

N 

N 

N 


53,2 - 43,1 

106141.400 

42.1 

5.29 

23 

N 

N 

N 

N 

N 


52,3 - 42,2 

106541.680 

27.2 

7.03 

23 

N 

N 

N 

Y 

N 


5i,4 - 4i,3 

109753.503 

18.8 

8.78 

22 

N 

N 

N 

N 

N 


61.5 - 5 i ,4 

131617.902 

25.1 

15.6 

19 

Y 

— 

N 

N 

N 


7i.7 - 61,6 

142701.325 

30.4 

20.2 

17 

Y 

— 

Y 

Y 

N 


7o,7 - 60,6 

146871.475 

28.3 

22.5 

17 

B 

— 

B 

B 

B 

CH 3 OCH 3 

72,6 - 62,5 

148223.143 

40.4 

21.2 

17 

B 

— 

N 

N 

B 

HCNH+ 

76,1 - 66,0 

148555.852 

135.7 

6.18 

17 

N 

— 

N 

N 

N 


76,2 - 66,1 

148555.852 

135.7 

6.18 

17 

N 

— 

N 

N 

N 


75,3 - 65,2 

148566.822 

103.0 

11.4 

17 

N 

— 

N 

N 

N 


75,2 - 65,1 

148566.823 

103.0 

11.4 

17 

N 

— 

N 

N 

N 


74,4 - 64,3 

148598.970 

76.2 

15.7 

17 

N 

— 

N 

N 

N 


74,3 - 64,2 

148599.354 

76.2 

15.7 

17 

N 

— 

N 

N 

N 


73,5 - 63,4 

148667.301 

55.3 

19.1 

17 

Y 

— 

N 

N 

N 


73,4 - 63,3 

148709.018 

55.4 

19.1 

17 

Y 

— 

Y 

N 

N 


76,2 - 66,1 

149792.574 

40.6 

21.9 

16 

N 

— 

N 

N 

Y* 


7i,6 - 61,5 

153432.176 

32.5 

25.1 

16 

Y 

— 

N 

N 

N 


132,11 - 13i,12 

155894.300 

105.9 

1.26 

16 

N 

— 

N 

N 

N 


122,10 - 12 i,ii 

155934.098 

92.4 

1.23 

16 

N 

— 

N 

N 

N 


ll2,9 - 111 ,10 

157072.457 

79.9 

1.22 

16 

N 

— 

N 

N 

N 


142,12 - 14i,13 

157115.035 

120.5 

1.32 

16 

N 

— 

N 

N 

N 


102,8 ^ 10l,9 

159127.569 

68.4 

1.21 

15 

— 

— 

N 

N 

N 


152,13 - 15i,14 

159739.080 

136.2 

1.39 

15 

— 

— 

N 

N 

N 


92,7 - 9i,8 

161899.774 

58.1 

1.22 

15 

— 

— 

N 

N 

N 


81,8 - 7 i ,7 

162958.657 

38.2 

30.5 

15 

— 

— 

N 

Y* 

Y 


82,6 - 81,7 

165176.756 

48.8 

1.22 

15 

— 

— 

N 

N 

N 


80,8 - 7 o ,7 

167320.697 

36.4 

33.5 

15 

— 

— 

N 

Y* 

Y 


72,5 - 7i,6 

168741.408 

40.6 

1.23 

15 

— 

— 

N 

N 

N 


82,7 - 72,6 

169299.154 

48.5 

32.6 

15 

— 

— 

N 

N 

N 


86,2 - 76,1 

169790.683 

143.9 

15.3 

14 

— 

— 

N 

N 

N 


86,3 - 76,2 

169790.683 

143.9 

15.3 

14 

— 

— 

N 

N 

N 


85,4 - 75,3 

169810.709 

111.1 

21.4 

14 

— 

— 

N 

N 

N 


85,3 - 75,2 

169810.715 

111.1 

21.4 

14 

— 

— 

N 

N 

N 


84,5 - 74,4 

169861.469 

84.3 

26.3 

14 

— 

— 

N 

N 

N 


84,4 - 74,3 

169862.523 

84.3 

26.3 

14 

— 

— 

N 

N 

N 


83,6 - 73,5 

169955.835 

63.5 

30.2 

14 

— 

— 

N 

N 

N 


83,5 - 73,4 

170039.076 

63.5 

30.3 

14 

— 

— 

N 

N 

Y 


82,6 - 72,5 

171620.760 

48.8 

33.9 

14 

— 

— 

N 

N 

N 


62,4 - 61,5 

172381.012 

33.4 

1.24 

14 

— 

— 

N 

N 

N 


lOl.io - 9 i ,9 

203335.761 

56.8 

60.3 

12 

Y 

— 

N 

N 

Y* 


lOo.lO - 00,9 

207679.189 

55.3 

64.7 

12 

Y 

— 

N 

Y* 

N 


102,9 - 02,8 

211328.960 

67.8 

65.6 

12 

Y 

— 

Y 

Y 

Yb 


105,6 - 95,6 

212323.555 

130.5 

52.0 

12 

N 

— 

N 

N 

N 


105,5 - 95,4 

212323.555 

130.5 

52.0 

12 

N 

— 

N 

N 

N 


104,7 - 94,6 

212428.020 

103.7 

58.4 

12 

N 

— 

N 

N 

B 

NH 2 CHO 

104,6 - 94,5 

212433.449 

103.7 

58.4 

12 

N 

— 

N 

N 

B 

NH 2 CHO 

103,8 ^ 93,7 

212572.837 

82.9 

63.3 

12 

Y 

— 

Y 

N 

N 
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Table B1 - continued 


Transition 

u 

(MHz) 

Eu 

(K) 

0 

1 

(") 

OMC-2 

CepE 

SVS13A 

IRAS4A 

116293 

Blends 

103,7 - 

03,6 

212832.307 

82.9 

63.6 

12 

N 

— 

N 

N 

N 


102,8 - 

02,7 

215687.009 

68.4 

69.8 

11 

Y 

— 

N 

N 

Y 


10l,9 - 

01,8 

218459.213 

60.8 

74.7 

11 

B 

— 

N 

N 

N 

CH 3 OH 

lll,ll - 

10 i,io 

223452.512 

67.5 

80.5 

11 

Y 

— 

N 

N 

N 


llo,ii - 

10 o,io 

227605.658 

66.2 

85.5 

11 

Y 

— 

N 

N 

N 


112,10 - 

102,9 

232273.646 

78.9 

88.2 

11 

N 

— 

N 

N 

Y 


ll5,7 - 

105,6 

233594.501 

141.7 

73.6 

11 

N 

— 

B 

N 

N 


ll5,6 - 

105,5 

233594.501 

141.7 

73.6 

11 

N 

— 

B 

N 

N 


ll4,8 - 

104,7 

233734.724 

114.9 

80.7 

11 

N 

— 

Y 

N 

N 


ll4,7 - 

104,6 

233745.613 

114.9 

80.7 

11 

N 

— 

Y 

N 

N 


ll3,9 - 

103,8 

233896.577 

94.1 

86.2 

11 

Y 

— 

Y 

N 

Y 


ll3,8 - 

103,7 

234315.498 

94.2 

86.7 

10 

— 

— 

N 

N 

N 


ll2,9 ^ 

102,8 

237896.684 

79.9 

94.8 

10 

Y 

— 

Y 

N 

B 

? 

111,10 - 

10l,9 

239951.800 

72.3 

99.6 

10 

Y 

— 

Y 

N 

N 


12i,12 - 

111,11 

243521.044 

79.2 

105 

10 

N 

— 

B 

N 

B 

CH 2 D 0 H 

12o,12 - 

llo,ii 

247390.719 

78.1 

110 

10 

N 

— 

N 

N 

N 


122,11 - 

112,10 

253165.793 

91.1 

115 

10 

Y 

— 

N 

N 

N 


124,9 - 

114,8 

255058.533 

127.2 

108 

10 

Y* 

— 

Y* 

N 

N 


124,8 - 

114,7 

255078.912 

127.2 

108 

10 

N 

— 

N 

N 

N 


123,10 - 

113,9 

255225.651 

106.4 

114 

10 

Y 

— 

Y 

N 

N 


123,9 - 

113,8 

255871.830 

106.4 

115 

10 

Y 

— 

N 

N 

N 


122,10 - 

ll2,9 

260189.090 

92.4 

125 

9 

B 

— 

B 

N 

B 

H 2 C 20 

12i,ii - 

111,10 

261327.450 

84.9 

129 

9 

N 

— 

Y 

N 

N 


13i,i3 - 

12i,12 

263542.236 

91.8 

133 

9 

Y 

— 

Y 

N 

N 



^Y: Detected above Tmb = 3cr. Y*: Weakly detected (S/N ~ 2 — 3; see Sect. [TTll . N: undetected. B: possibly detected but blended, 
not observed. 

^Detected but with an anomalously high flux (maybe blended); removed from analysis. 





Table B2. HNCO transitions searched for in this study and 3(T detections' 


Transition 

V 

(MHz) 

Eu 

(K) 

^ul 

(10-^ s-1) 

(") 

OMC-2 

CepE 

SVS13A 

IRAS4A 

116293 

L1157 

L1527 

B1 

LI 544 

TMC-1 

4i,4 - 

3i,3 

87597.330 

53.8 

0.80 

28 

N 

N 

Y 

Yb 

N 

N 

N 

N 

N 

— 

4o,4 - 

3o,3 

87925.237 

10.5 

0.88 

28 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

— 

4i,3 - 

3i,2 

88239.020 

53.9 

0.82 

28 

Y* 

N 

■y"* 

Yb 

N 

N 

N 

N 

N 

— 

5i.5 - 

4i,4 

109495.996 

59.0 

1.7 

22 

N 

N 

Y 

N 

Y 

W 

N 

N 

N 

N 

5o,5 - 

4o,4 

109905.749 

15.8 

1.8 

22 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

5i,4 - 

4i,3 

110298.089 

59.2 

1.7 

22 

Y* 

N 

Y 

N 

Y 

N 

N 

N 

N 

N 

61,6 - 

5i,5 

131394.230 

65.3 

2.9 

19 

N 

N 

Y 

N 

Y 

N 

N 

N 

— 

N 

60.6 - 

5o,5 

131885.734 

22.2 

3.1 

19 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

— 

Y 

61.5 - 

5i,4 

132356.701 

65.5 

3.0 

19 

N 

N 

N 

N 

Y 

Y 

N 

N 

— 

N 

7i,7 - 

61,6 

153291.935 

72.7 

4.7 

16 

N 

— 

Y 

N 

Y 

N 

N 

N 

— 

N 

7o,7 - 

60,6 

153865.086 

29.5 

4.9 

16 

Y 

— 

Y 

Y 

Y 

Y 

Y 

Y 

— 

Y 

7i,6 - 

6 i,s 

154414.765 

72.9 

4.8 

16 

N 

— 

Y 

N 

Y 

N 

N 

N 

— 

N 

10i,io 

- 9i,9 

218981.009 

101.1 

14.2 

11 

N 

N 

Y 

N 

Y 

N 

N 

N 

— 

— 

lOo.io 

- 9o,9 

219798.274 

58.0 

14.7 

11 

Y 

Y 

Y 

Y 

Y 

N 

N 

Y 

— 

— 

101.9 - 

- 9i,8 

220584.751 

101.5 

14.5 

11 

Y 

N 

Y 

Y 

Y 

N 

N 

N 

— 

— 

111,11 - 

10i,io 

240875.727 

112.6 

19.0 

10 

N 

— 

Y 

Y 

Y 

N 

N 

N 

— 

— 

llo.ll - 

10o,io 

241774.032 

69.6 

19.6 

10 

B 

— 

B 

B 

B 

N 

N 

N 

— 

— 

111,10 - 

10l,9 

242639.704 

113.1 

19.5 

10 

N 

— 

Y 

N 

N 

N 

N 

— 

— 


12i,i2 - 

111,11 

262769.477 

125.3 

24.8 

9 

N 

N 

Y 

Y 

yc 

N 

N 

N 

— 

— 

12o,12 - 

llo,ll 

263748.625 

82.3 

25.6 

9 

Y 

Y 

Y 

Y 

Y 

N 

N 

N 

— 

— 

12i,ii - 

111,10 

264693.655 

125.9 

25.4 

9 

N 

N 

Y 

Y 

Y 

N 

N 

N 

— 

— 


^Y: Detected above Tmh = 3cr. Y*: Weakly detected (S/N 2 — 3; see Sect. [XTl . N: undetected. B: detected but blended. —: not observed. 
'^Detected but with an anomalously high flux (maybe blended): removed from analysis. 

^Blended with an unidentified feature: removed from analysis. 

Blends: CH 3 OH at 241.774 GHz 
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Table B3. L1544: Gaussian fits to the detected HNCO lines 


Transition u RMS Tpeak Msr f 

(MHz) (mK) (mK) (km s“^) (km s“^) (mK km s“^) 

40.4 - 3o,3 87925.237 11.9 459 (2) 7.2 (1) 0.7 (2) 342 (9) 

50.5 - 4o,4 109905.996 4.8 601 (2) 7.6 (1) 0.7 (1) 448 (5) 


Table B4. TMC-1: Gaussian fits to the detected HNGO lines 


Transition i/ RMS Tpeak Msr f 

(MHz) (mK) (mK) (km s“^) (km s“^) (mK km s“^) 


5o,5 - 4o,4 

109905.996 

6.8 

157 (8) 

5.8 (1) 

1.2 (1) 

203 (11) 

6o,6 - 5o,5 

131885.734 

6.7 

94 (7) 

5.8 (1) 

0.5 (3) 

50 (5) 

7o,7 - 00,6 

153865.086 

3.2 

28 (3) 

5.9 (1) 

0.5 (4) 

15 (2) 


Table B5. Bl: Gaussian fits to the detected HNGO lines 


Transition 

V 

(MHz) 

RMS 

(mK) 

rn 

peak 

(mK) 

(km s“^) 

AH 

(km s“^) 

/ T^hdV 
(mK km s~^) 

4o,4 - 3o,3 

87925.237 

3.0 

530 (12) 

0.7 (1) 

1.4(1) 

765 (5) 

5o,5 - 4o,4 

109905.996 

22.2 

480 (6) 

6.6 (1) 

1.3 (1) 

662 (32) 

00,6 - 5o,5 

131885.734 

6.3 

345 (7) 

0.7(1) 

1.4(1) 

521 (8) 

7o,7 - 00,6 

153865.086 

7.4 

224 (6) 

6.6 (1) 

1.4(1) 

326 (10) 


Table B6. L1527: Gaussian fits to the detected HNGO lines 


Transition 

V 

(MHz) 

RMS 

(mK) 

^peak 

(mK) 

Vl.r 

(km s“^) 

AK 

(km s~^) 

/ T^hdV 
(mK km s“^) 

4o,4 - 3o,3 

87925.237 

2.3 

145 (2) 

5.9 (1) 

1.3 (1) 

198 (3) 

5o,5 - 4o,4 

109905.996 

7.9 

135 (7) 

5.9 (1) 

1.2 (1) 

175 (11) 

00,6 - 5o,5 

131885.734 

6.7 

138 (1) 

5.9 (1) 

0.8 (1) 

115 (8) 

7o,7 - 00,6 

153865.086 

8.0 

03 (1) 

5.8 (1) 

0.7 (1) 

47 (0) 


Table B7. L1157mm: Gaussian fits to the detected HNCO lines 


Transition 

V 

(MHz) 

RMS 

(mK) 

^peak 

(mK) 

(km s“^) 

AK 

(km s~^) 

/ T^hdV 
(mK km s“^) 

4o,4 - 3o,3 

87925.237 

3.3 

113 (7) 

2.0 (1) 

1.0 (1) 

198 (7) 

5o,5 - 4o,4 

109905.996 

7.8 

142 (5) 

2.5 (1) 

1.2 (1) 

177 (13) 

00,6 - 5o,5 

131885.734 

5.5 

81 (3) 

2.0 (1) 

1.5 (1) 

126 (11) 

7o,7 - 00,6 

153865.086 

5.3 

71 (1) 

2.0 (1) 

1.1 (1) 

81 (8) 
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Table B8. IRAS 4A: Gaussian fits to the detected NH 2 CHO and HNCO lines 


Transition 

V 

RMS 

^peak 


AH 



(MHz) 

(mK) 

(mK) 

(km s~^) 

(km s~l) 

(mK km s~^) 

NH2CH0 

5o,5 - 4o,4 

105464.219 

3.3 

9.8 (1) 

6.7 (4) 

4.0 (13) 

41 (11) 

52,3 - 42,2 

106541.680 

2.2 

7.7 (17) 

7.6 (3) 

2.8 (9) 

23 (5) 

7i,7 - 61,6 

142701.325 

5.0 

15 (4) 

8.1 (3) 

2.6 (5) 

40 (8) 

00 

00 

1 

162958.657 

38.2 

14 (6) 

8.5 (4) 

2.9 (7) 

44 (9) 

00 

0 

00 

t 

0 

167320.697 

36.4 

13 (4) 

7.2 (4) 

3.2 (7) 

43 (10) 

10o,io - 9 o,9'^ 

207679.189 

55.3 

19 (10) 

6.6 (5) 

3.3 (10) 

66 (19) 

102,9 - 92,8 

211328.960 

7.1 

31 (3) 

6.7 (2) 

2.4 (5) 

82 (15) 

HNCO 

4o,4 - 3o,3 

87925.237 

3.1 

195 (16) 

7.2 (1) 

2.2 (1) 

458 (6) 

5o,5 - 4o,4 

109905.749 

8.0 

198 (15) 

7.1 (1) 

2.3 (1) 

495 (16) 

60,6 - 5o,5 

131885.734 

5.9 

203 (12) 

7.0 (1) 

2.5 (1) 

545 (11) 

7o,7 - 60,6 

153865.086 

6.8 

168 (11) 

7.1 (1) 

2.2 (1) 

395 (12) 

10o,10 - 9 o,9 

219798.274 

9.0 

88 (8) 

6.8 (1) 

3.3 (3) 

307 (19) 

10l,9 - 9 i,8 

220584.751 

6.8 

37 (6) 

6.5 (3) 

2.8 (7) 

111 (22) 

111,11 - 10i,io 

240875.727 

7.2 

23 (2) 

6.6 (2) 

1.7 (6) 

42 (13) 

12i,i2 - 111,11 

262769.477 

8.9 

44 (4) 

6.6 (2) 

3.8 (6) 

177 (21) 

12o, 12 - llo,ll 

263748.625 

12.9 

54 (4) 

6.5 (8) 

3.9 (29) 

230 (130) 

12i,ii - 111,10 

264693.655 

9.1 

26 (5) 

6.3 (4) 

3.7 (9) 

103 (21) 


^ Transition weakly detected (see Table [BT} but included in the analysis for completeness. 
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Table B9. 116293: Gaussian fits to the detected NH2CHO and HNCO lines (intensity in units) 


Transition 

u 

RMS 

'T 

peak 

Vi„ 

AH 

/ T,dV 


(MHz) 

(mK) 

(mK) 

(km s~^) 

(km s“^) 

(mK km s“^) 

NH 2 CH 0 

4o,4 - 3o,3 

84542.330 

5.3 

17 (6) 

3.2 (12) 

6.7 (20) 

120 (40) 

4i,3 - 3i,2 

87848.873 

2.6 

12 (4) 

2.6 (5) 

3.8 (11) 

50 (13) 

5i,5 - 4i_4 

102064.267 

3.3 

16 (5) 

2.0 (3) 

2.4 (9) 

42 (12) 

5o,5 - 4o,4 

105464.219 

5.3 

20 (7) 

5.1 (4) 

3.0 (10) 

64 (19) 

76,2 - 66 , 1 '" 

149792.574 

10.0 

29 (10) 

2.6 (4) 

3.0 (7) 

93 (22) 

81,8 - 7 i ,7 

162958.657 

10.1 

26 ( 10 ) 

2.8 (4) 

2.5 (13) 

69 (27) 

00 

0 

00 

1 

0 

167320.697 

10.6 

55 (11) 

2.0 ( 2 ) 

2.3 (5) 

136 (24) 

83,5 - 73,4 

170039.076 

15.3 

50 (80) 

2.7 (17) 

2.4 (46) 

120 (170) 

lOi.io - 9 i, 9 "' 

203335.761 

6.8 

19 (7) 

0.5 (10) 

6.7 (27) 

135 (43) 

102,8 - 92,7 

215687.009 

5.5 

30 (10) 

1.0 (13) 

5.5 (34) 

175 (90) 

112,10 - 102,9 

232273.646 

4.6 

23 (18) 

1.5 (19) 

5.0 (39) 

120 (90) 

ll2,9 - 102,8 

237896.684 

8.2 

45 (9) 

2.8 (3) 

3.1 (15) 

148 (42) 

HNCO 

4o,4 - 3o,3 

87925.237 

3.2 

162 (3) 

4.0 (1) 

3.9 (1) 

671 (11) 

5i,5 - 4i,4 

109495.996 

5.4 

38 ( 6 ) 

2.0 (3) 

4.5 (11) 

182 (35) 

5o,5 - 4o,4 

109905.749 

5.9 

254 ( 6 ) 

3.9 (1) 

3.2 (1) 

853 (18) 

5i,4 - 4i,3 

110298.089 

5.4 

16 (5) 

5.4 ( 6 ) 

5.4 (18) 

90 (23) 

61,6 - 5i,s 

131394.230 

5.4 

39 (17) 

2.5 ( 8 ) 

5.2 (27) 

217 ( 88 ) 

60,6 - 5 o,s 

131885.734 

7.3 

263 ( 8 ) 

3.9 (1) 

3.0 (1) 

825 (18) 

61,5 ^ 5 i ,4 

132356.701 

6.9 

39 ( 8 ) 

3.6 (5) 

5.3 (13) 

218 (48) 

7i,7 - 61,6 

153291.935 

9.5 

73 (27) 

5.1 (5) 

3.0 (16) 

230 (90) 

7o,7 - 60,6 

153865.086 

12.3 

224 (18) 

3.9 (1) 

4.0 (3) 

950 (50) 

7 i ,6 - 61.S 

154414.765 

11.1 

57 (28) 

2.7 (7) 

3.4 (21) 

200 (90) 

10 i,io - 9 i ,9 

218981.009 

6.2 

103 ( 8 ) 

2.8 ( 2 ) 

6.0 (5) 

664 (47) 

10o,10 - 9o,9 

219798.274 

4.4 

260 ( 12 ) 

3.3 (1) 

5.7 (3) 

1580 (60) 

10l,9 - 9 i,8 

220584.751 

6.4 

85 (10) 

2.9 (3) 

6.2 (7) 

560 (50) 

111,11 - 10 i,io 

240875.727 

14.5 

132 (16) 

2.8 (3) 

6.6 (9) 

930 (100) 

12o, 12 - llo,ll 

263748.625 

7.9 

270 (70) 

3.7 ( 6 ) 

5.6 (15) 

1640 (360) 

12141 - 111,10 

264693.655 

6.9 

119 (11) 

3.4 (2) 

6.6 ( 6 ) 

840 (60) 


Transition weakly detected (see Table [bh but included in the analysis for completeness. 
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Table BIO. SVS13A: Gaussian fits to the detected NH2CHO and HNCO lines 


Transition 

V 

RMS 

^peak 

V)sr 

AH 



(MHz) 

(mK) 

(mK) 

(km s“l) 

(km s“l) 

(mK km s~^) 

NH2CH0 

5o,5 - 4o,4 

105464.219 

3.9 

16 (1) 

8.1 (2) 

1.8 (5) 

14 (1) 

7i,7 - 61,6 

142701.325 

6.4 

28 (3) 

8.6 (2) 

2.2 (5) 

22 (5) 

73,4 - 63,3 

148709.018 

7.1 

26 (2) 

7.1 (3) 

2.2 (7) 

22 (3) 

102,9 - 92,8 

211328.960 

7.3 

62 (5) 

7.5 (5) 

4.4 (14) 

31 (2) 

103,8 - 93,7 

212572.837 

7.2 

43 (7) 

8.4 (2) 

2.7 (6) 

29 (2) 

114,8 - 104,7 

233734.724 

6.6 

26 (2) 

6.9 (4) 

3.0 (9) 

23 (5) 

114,7 - 104,6 

233745.613 

7.5 

26 (2) 

8.2 (3) 

1.5 (7) 

26 (1) 

113,9 ^ 103,8 

233896.577 

6.6 

49 (14) 

8.2 (4) 

2.9 (9) 

20 (3) 

ll2,9 - 102,8 

237896.684 

6.9 

48 (6) 

8.1 (3) 

4.3 (8) 

28 (4) 

111,10 - 10l,9 

239951.800 

8.1 

66 (7) 

7.6 (4) 

3.4 (9) 

23 (2) 

123,10 - 113,9 

255225.651 

5.6 

45 (2) 

8.8 (12) 

3.3 (36) 

24 (4) 

12l,ll - 111,10 

261327.450 

8.7 

26 (2) 

8.4 (2) 

4.6 (6) 

38 (5) 

13i,13 - 12i,12 

263542.236 

7.6 

54 (5) 

7.7 (6) 

4.4 (15) 

32 (5) 

HNCO 

4i,4 - 3i,3 

87597.330 

3.5 

11 (3) 

6.4 (4) 

3.3 (8) 

40 (9) 

4o,4 - 3o,3 

87925.237 

9.1 

46 (4) 

8.6 ( 2 ) 

3.1 (6) 

155 (23) 

4i,3 - 31,2" 

88239.020 

3.0 

13 (3) 

7.3 (4) 

3.3 (20) 

45 (17) 

5i,5 - 4i,4 

109495.996 

3.9 

16 (1) 

7.1 (2) 

2.0 (5) 

35 (8) 

5o,5 - 4o,4* 

109905.749 

5.4 

76 (5) 

8.5 (1) 

1.2 ( 1 ) 

100 ( 8 ) 

5i,4 - 4i,3 

110298.089 

6.2 

28 (4) 

8.2 (3) 

4.7 (8) 

139 (19) 

61,6 - 5i,s 

131394.230 

6.6 

19 (3) 

8.2 (4) 

3.6 (8) 

72 (15) 

60,6 - 5o,5* 

131885.734 

6.0 

77 (10) 

8.5 (1) 

1.6 (3) 

129 (15) 

7i,7 - 61,6 

153291.935 

5.7 

25 (2) 

8.5 (5) 

5.1 (13) 

140 (30) 

7o,7 - 60,6* 

153865.086 

6.1 

70 (6) 

7.5 (1) 

3.4 (5) 

258 (19) 

7i,6 - 61,5 

154414.765 

6.0 

32 (3) 

8.3 (6) 

4.4 (15) 

149 (40) 

10i,io - 9i,9 

218981.009 

6.8 

46 (3) 

8.2 ( 2 ) 

4.3 (4) 

212 (15) 

10 o ,10 - 9 o ,9 

219798.274 

7.1 

89 (4) 

8.4 (1) 

3.2 (2) 

309 (16) 

10l,9 - 9i,8 

220584.751 

6.0 

35 (3) 

8.5 (4) 

3.7 (9) 

136 (27) 

111,11 - 10i,io 

240875.727 

6.7 

46 (4) 

8.2 (5) 

4.9 (17) 

236 (66) 

111,10 - 10l,9 

242639.704 

8.7 

42 (5) 

8.1 (3) 

4.6 (6) 

206 (24) 

12i,i2 - 111,11 

262769.477 

9.9 

68 ( 4 ) 

8.5 (2) 

5.3 (5) 

380 (32) 

12o,12 - llo,ll 

263748.625 

9.3 

60 (5) 

8.2 (4) 

3.4 (9) 

213 (48) 

12 i,ii - 111,10 

264693.655 

8.9 

39 (4) 

8.5 (3) 

4.4 ( 6 ) 

183 (23) 


* Transition affected by emission at OFF position: lower limit point in the rotational diagram. 
Transition weakly detected (see Table [bh but included in the analysis for completeness. 


Table Bll. Cep E: Gaussian fits to the detected NH2CHO and HNCO lines 


Transition 

V 

RMS 

^peak 

Visr 

AH 



(MHz) 

(mK) 

(mK) 

(km s“^) 

(km s”^) 

(mK km s“^) 

NH2CH0 

4o,4 - 3o,3 

84542.330 

1.9 

7.7 (1) 

- 11.5 (4) 

2.8 (10) 

23 (7) 

42,2 - 32,1 

85093.272 

1.3 

4.3 (1) 

-11.8 (3) 

3.5 (7) 

16 (3) 

41,3-31,2" 

87848.873 

1.2 

6.3 (1) 

-10.5 (2) 

2.6 (5) 

17 (3) 

5i,5 - 4i,4 

102064.267 

1.6 

7.6 (1) 

-10.6 (2) 

1.7 (5) 

13 (3) 

5o,5 - 4o,4 

105464.219 

2.1 

4.1 (1) 

-10.6 (6) 

2.9 (11) 

13 (4) 

HNCO 

4o,4 - 3o,3 

87925.237 

1.4 

90 (4) 

-11.1 (1) 

1.9 (1) 

179 (3) 

5o,5 - 4o,4 

109905.749 

2.8 

104 (12) 

-11.1 (1) 

2.4(1) 

262 (6) 

60,6 - 5o,s 

131885.734 

7.8 

128 (9) 

-11.1 (1) 

2.2 (3) 

299 (23) 

10o,10 - 9 o,9 

219798.274 

7.5 

45 (10) 

-10.1 (73) 

8.0 (80) 

380 (65) 

12o, 12 - llo,ll 

263748.625 

6.0 

38 (8) 

-10.2 (5) 

5.6 (14) 

226 (45) 


Transition weakly detected (see Table [511 but included in the analysis for completeness. 
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Table B12. OMC-2 FIR 4: Gaussian fits to the detected NH2CHO and HNCO lines 


Transition 

V 

RMS 

^peak 


AH 



(MHz) 

(mK) 

(mK) 

(km s~^) 

(km s~l) 

(mK km s~^) 

NH2CH0 

4o,4 - 3o,3 

84542.330 

2.8 

15 (2) 

10.8 (2) 

2.2 (4) 

35 (5) 

4i,3 - 3i,2 

87848.873 

2.8 

9 (2) 

11.0 (4) 

3.6 (6) 

37 (7) 

5i,5 - 

102064.267 

3.5 

15 (1) 

11.0 (2) 

2.9 (5) 

46 (6) 

5o,5 - 4o,4 

105464.219 

4.7 

16 (1) 

11.5 (4) 

3.7 (10) 

64 (15) 

52,4 - 42,3 

105972.599 

4.6 

15 (5) 

11.3 (2) 

1.7 (5) 

28 (8) 

61,5 - 5i,4 

131617.902 

5.6 

29 (3) 

11.2 (2) 

3.8 (6) 

115 (13) 

7i,7 - 61,6 

142701.325 

6.4 

28 (3) 

11.5 (2) 

3.1 (7) 

91 (14) 

73,4 - 63,3 

148709.018 

5.6 

26 (2) 

11.5 (2) 

1.8 (5) 

51 (10) 

7i,6 - 61,5 

153432.176 

8.2 

40 (3) 

11.8 (3) 

2.7 (10) 

114 (31) 

10i,io - 9 i,9 

203335.761 

10.8 

52 (4) 

11.3 (2) 

2.8 (5) 

156 (20) 

10o,10 - 9o,9 

207679.189 

8.8 

49 (3) 

11.6 (1) 

1.8 (4) 

97 (14) 

103,8 ^ 93,7 

212572.837 

12.7 

43 (7) 

11.7 (3) 

2.9 (6) 

131 (23) 

111,11 - 10i,io 

223452.512 

12.6 

44 (6) 

11.7 (5) 

3.0 (12) 

142 (48) 

llo,ii - 10o,io 

227605.658 

14.5 

61 (10) 

11.9 (3) 

3.7 (7) 

237 (33) 

113,9 - 103,8 

233896.577 

16.9 

49 (14) 

11.6 (3) 

2.9 (10) 

151 (36) 

ll2,9 - 102,8 

237896.684 

10.1 

48 (6) 

11.4 (2) 

3.7 (6) 

186 (23) 

111,10 - 10l,9 

239951.800 

10.5 

66 (7) 

11.5 (1) 

2.5 (4) 

175 (21) 

122,11 - 112,10 

253165.793 

12.4 

44 (9) 

11.2 (3) 

3.2 (7) 

152 (27) 

124.9 - 114,8'" 

255058.533 

12.5 

39 (12) 

11.6 (2) 

1.9 (4) 

80 (17) 

123,10 - 113,9 

255225.651 

11.4 

45 (2) 

11.5 (4) 

1.6 (9) 

74 (4) 

123,9 - 113,8 

255871.830 

11.8 

41 (7) 

11.8 (2) 

2.7 (5) 

119 (19) 

HNCO 

4o,4 - 3o,3 

87925.237 

2.8 

128 (13) 

11.2 (1) 

2.5 (1) 

344 (6) 

4i,3 - 3i,2'" 

88239.020 

2.3 

9 (2) 

12.7 (3) 

3.1 (6) 

31 (5) 

5o,5 - 4o,4 

109905.749 

6.9 

224 (22) 

11.4 (1) 

2.1 (1) 

512 (14) 

5i,4 -4i,3'" 

110298.089 

6.2 

18 (5) 

12.8 (2) 

1.4 (5) 

26 (7) 

60,6 - 5o,s 

131885.734 

4.7 

267 (27) 

11.2 (1) 

2.7 (1) 

775 (11) 

7o,7 - 60,6 

153865.086 

8.8 

304 (38) 

11.2 (1) 

3.0 (4) 

962 (89) 

lOo.lO - 9 o ,9 

219798.274 

9.2 

343 (28) 

11.4 (1) 

2.5 (1) 

918 (17) 

10l,9 - 9 i,8 

220584.751 

11.0 

36 (1) 

11.7 (4) 

2.1 (13) 

79 (35) 

12o,12 - llo.ll 

263748.625 

5.1 

202 (9) 

11.9 (3) 

3.0 (5) 

648 (24) 


Transition weakly detected (see Table [bH but included in the analysis for completeness. 
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22 Lopez-Sepulcre et al. 
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Figure Cl. HNCO observed spectral lines (black) in L1544, TMC-1, Bl, L1527, and L1157mm, and the spectra predicted by best fit 
LTE model (red). 
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Figure C2. Sample of HNCO (/e/i) and NH 2 CHO {right) observed spectral lines (black) in IRAS 4A and SVS13A (compact solution), 
and the spectra predicted by best fit LTE model (red). 
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Figure C3. Sample of HNCO (left) and NH2CHO (right) observed spectral lines (black) in Cep E and OMC-2 FIR 4 (extended 
solutions), and the spectra predicted by best fit LTE model (red). 
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Figure C4. Rotational diagrams for L1544, TMC-1, L1527, and L1157mm IRAS 4A (middle), and 116293 (right). Data points 

are depicted in black. The red lines correspond to the best fit to the data points. The dashed vertical lines in the middle and right panels 
indicate the upper-level energy (35 K) at which the division of the 2-component fitting was made. 
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Figure C5. Rotational diagrams for SVS13A (left), OMC-2 FIR 4 (middle), and Cep E (right). Data points are depicted in black. The 
red lines correspond to the best fit to the data points. . 





































